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1. Introduction and motivation 

Despite the huge efforts of the last years, up to now, we still do not have a clear understanding 
about the reasons why quarks and gluons remain confined in colour singlets. 

The infrared properties of the gluon and ghost propagators in Landau gauge have been related 
with possible gluon confinement mechanisms. In this sense, the computation of the propagators 
can help in the understanding of the QCD confinement mechanism [1]. 

1.1 Gluon confinement in Landau gauge 

The so-called Kugo-Ojima scenario (KO) [§] and the Gribov-Zwanziger horizon condition 
(GZ) [||, 0] impose restrictions on the infrared behaviour of the Landau gauge gluon and ghost 
propagators, 

DfM = 5"' [d,. - ^-f^^ Dig') , (LI) 
G"'{q) = -6"'G{q^), (1.2) 

The GZ mechanism requires D{0) =0 (which implies maximal violation of reflection positiv- 
ity) and an enhanced ghost propagator, relative to the perturbative function. The KO confinement 
mechanism demands 1 / {q^G{q^)) = in the limit q ^ Q. From the point of view of the KO and 
GZ confinement mechanisms, the requirements on D(0) and G(0) are necessary conditions and its 
violation immediatly rules out these scenarios. 

The infrared behaviour of these propagators can not be studied using perturbative methods. 
One should rely on non-perturbative methods. One possibility is the use of Dyson-Schwinger 
equations, which are an infinite tower of coupled nonlinear equations for the QCD Green's func- 
tions. Although they can provide an analytical solution in the infrared, the computation of the 
gluon and ghost propagators within DSE requires defining a truncation scheme and parameterizing 
some of the QCD vertices. 

For the Landau gauge, recent solutions can be classified in two categories: solutions with 
a finite D(0), which are not compatible with the KO or GZ confining mechanism (see [Q] and 
references therein) and solutions which do not rule out the two confining mechanisms (see [§] for 
a review). For the later class of solutions, they predict a pure power law behaviour in the infrared, 

Z{q^) = q^D{q^)=C0{q'f\ (1.3) 
F{q^) = q^G{q^)=C0'{qY"\ (1-4) 

with K = 0.595, which, for the zero momentum, implies a null (infinite) gluon (ghost) propagator. 

Lattice QCD simulations provide another method of study these quantities non perturbatively, 
although one should care about finite volume and finite lattice spacing effects. Despite efforts 
in large simulations in recents years (see, for example, [0, ||]), lattice results show a finite non- 
vanishing zero momentum gluon propagator, with no sign for a turnover in the infrared region. 

In this paper we present an update on our study of Cucchieri-Mendes bounds [Q] in SU(3) 
lattice gauge theory [^. 
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1.2 Cucchieri-Mendes bounds 

The Cucchieri-Mendes bounds provide upper and lower bounds for the zero momentum 
gluon propagator of lattice Yang-Mills theories in terms of the average value of the gluon field. In 
particular, they relate the gluon propagator at zero momentum D{0) with 

^(0) = ^7(]^IK(0)I' (1-5) 



where d is the number of space-time dimensions, and Nc the number of colors. In the above 
equation, (0) is the a color component of the gluon field at zero momentum, defined by 

^m(O) = ^I^mW (1-6) 

where A^(x) is the a color component of the gluon field in the real space. According to [|8|], the 
D(0) is related with M(0) by 

(M(0))2 < M < ^ (A^2 _ 1) ^^(0)2^ _ (17) 
In the last equation ( ) means Monte Carlo average over gauge configurations. The bounds in 



equation ( 1 .7 ) are a direct result of the Monte Carlo approach. The interest on these bounds comes 



from allowing a scaling analysis which can help understanding the finite volume behaviour of 



D{0): assuming that each of the terms in inequality (1.7) scales with the volume according to 
A/V", the simplest possibility and the one considered in [j8|], an a > 1 for (M(0)^) clearly indicates 
that D{0) —^0 as the infinite volume is approached. In this sense, this scaling analysis allows to 
investigate the behaviour of D{0) in the infinite volume limit. 

For the SU(2) Yang-Mills theory |j8|], the results show a D(0) =0 for the two dimensional 
theory, but a D{0) ^ for three and four dimensional formulations. Our analysis for SU(3) favors 
a D(0) = result. 



2. Cucchieri-Mendes bounds in SU(3) gauge theory 

We have studied the Cucchieri-Mendes bounds for SU(3) gauge theory with two values of the 
gauge coupling: j8 = 6.0 and j8 = 5.7. 

2.1 Scaling analysis for j3 = 6.0 

For j3 =6.0, we used the lattice setup described in table [T[ The differences with [^] being that 
we now use a new ensemble for 24^^, and we also use preliminary data for 80^ lattices. 

In figure % we show (M(0))^ D(0)/y and d{Nl - 1)(M(0)2), together with the fits to 
In table ^ we see the results of fitting the lattice data to A/V". The values for a strongly 
support a vanishing zero momentum gluon propagator. However, the fits reported in table ^ to 
CjW + ft)V^", support a D(0) = C 7^ 0. So, no conclusive answer can be made from the analysis 
of the two fits. 
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Figure 1: Cucchieri-Mendes bounds for the lattice data at j3 = 6.0. 



The fact that the a's in table ^ are in disagreement with the results for SU(2) can have 
several explanations. First of all, there can be differences between calculations done using different 
gauge groups, although recent studies seem to support the equivalence of the two theories ^ [10, n|]. 

Other possibilities for such a difference are the use of different lattice volumes and/or lattice 
spacing. Indeed, the physical volumes considered in the SU(2) study, up to [27 fm)^, are larger 



Z? l¥ W W 28^^ 3¥ 48^ 64^ 80^ 
L(fm) 1.63 2.03 2.44 2.84 3.25 4.88 6.50 8.13 
#conf 52 12 60* 56 126 104 120 18* 

Table 1: Lattice setup for j3 = 6.0 configurations. The lattice spacing is a = 0. 1016(25)fm. * stands for new 
data, comparing with . 





CO 


a 


Xv 


(M(0)) 


9.66(37) 


0.5265(27) 


0.69 




151±11 


1.0551(51) 


0.55 


NM-^)m^?) 


2999 ±228 


1.0522(54) 


0.71 



Table 2: Fits toA/V" using lattice data at j3 = 6.0. 





CO 


a 


C 


Xv 


{Mm' 


125 ±74 


1.148(8) 


28.2±8.9 


0.57 


m/v 


247 ±201 


1.17(10) 


46 ±13 


0.40 


NM-^){M{Of) 


5215 ±258 


1.18(11) 


1015 ±258 


0.60 



Table 3: Fits to C/V + (oV^" — lattice data at j3 = 6.0. 



Note that a recent direct comparison between SU(2) and SU(3) gluon propagators showed a measurable difference 
in the infrared region — [O]. 
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^4 84 184 36^ 

L(fm) 1.47 1.84 2.57 3.31 4.78 6.62 8.09 
#conf 56 149 149 149 132 100 29 



Table 4: Lattice setup for j3 = 5.7 configurations. The lattice spacing is a = 0. 1838(1 l)fm. 

than the ones used in our SU(3) study — up to (8fm)4. On the other hand, the lattice spacing used 
in the SU(2) study is more than twice our lattice spacing. 

2.2 Scaling analysis for j3 =5.7 

To avoid considering larger volumes, which is computationally demanding, we carried simu- 
lations with a lower p value, in the case j3 = 5.7 (a = 0.1838(ll)fm), keeping volumes similar to 
the ones considered previously — see table The new lattice setup is reported in table 

In figure ^ we can see the comparison between the propagators computed using different lattice 
spacings. 

Note that, for the two lowest volumes, we have also plotted some data for p = 6.2 configu- 
rations. Although the j3 = 5.7 and 6.0 data seem to be similar, there are some differences, in the 
infrared, with the propagators computed at j3 = 6.2. This deserves further investigations to clarify 
any possible effects due to finite lattice spacing. 

In what concerns the Cucchieri-Mendes bounds, we can see our results in figure |3[ Moreover, 
in table ^ we see the results of fitting^ the lattice data toA/V". 

As in the p = 6.0 case, the a's values strongly support a D(0)=0. However, the lattice data 
is also compatible with the functional form C/V + (OV^'^ — see table ^ The figures presented in 
this table do not allow a conclusive answer, as the fit to D(0)/V allows for D(0) = 0; however, the 
C figure for {M{0)f imphes a D(0) / 0. 





CO 


a 


Xv 


(M(0)) 


4.73(13) 


0.5267(25) 


1.22 


m/v 


34.0 ±1.7 


1.0504(45) 


0.80 




726±41 


1.0530(50) 


1.08 



Table 5: Fits to A /V" using lattice data at j3 = 5.7. 





ft) 


a 


C 


Xv 


(M(0))2 


22 ±8 


1.13(10) 


7.0±3.7 


1.40 


m/v 


32 ±13 


1.063(91) 


3 ±22 


1.07 


Nd{N^-l){M{Oy-) 


688 ±189 


1.117(90) 


212±135 


1.29 



Table 6: Fits to C/V + ©V"" — lattice data at j3 = 5.7. 



^In the fits to j8 = 5.7 data, to keep Xv <^ the 26^ lattice data had to be excluded. 
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Renormalized Gluon Propagator - |J. = 3 GeV 



Renormalized Gluon Propagator - H = 3 GeV 



,15 



18 -p = 5.7 On Axis 
Conical [3.31 fm] 

Si''- 6.0 On Axis 
Conical [3.25 fml 

48"' - p ^ 6.2 On Axis 
Conical [3.49 fm] 



q [GeV] 

(a) L ~ 33 fm 



o 26"* 


- p ^ 5.7 On Axis 




Conical [4.78 fm] 


° 48'' 


- p ^ 6.0 On Axis 




Conical [4.88 fm] 


" 64'' 


-p^6.2 0n Axis 




Conical [4.65 fm] 



C| IGeV] 

(b) L ~ 4.8/m 



Renormalized Gluon Propagator - |J. = 3 GeV 



Renormalized Gluon Propagator - H = 3 GeV 



4 r 
3 I- 
2 h 



q [GeVl 

(c) L ~ 6.5 fm 



36 -p ^5.7 On Axis 
Conical [6.62 fm] 

64"' - p ^ 6.0 On Axis 
Conical [6.50 fm] 



q IGeV] 

(d) L~8.1/m 



44 - [1 ^ 5.7 On Axis 
Conical r8 09 fm] 

80^ - P - 6.0 On Axis 
Conical [8.13 fml 



Figure 2: Comparing the gluon propagator computed with different lattice spacings at the same physical 
volume. 



3. Conclusions and future work 



In this article we presented a study of Cucchieri-Mendes bounds in SU(3) pure gauge theory, 
using two different values of lattice spacing. Fitting the data to a pure power law in the volume 
A/V", one finds always D{0) = 0. Note that, although this result is in disagreement with the SU(2) 
study [|8]], it is in agreement with our study presented in [[l3|], where we computed a fc ~ 0.53. This 
value for k implies a vanishing zero momentum gluon propagator. 

However, the use of more general ansatze for the dependence with the lattice volume do not 
allow to take definitive conclusions. This puzzle deserves a more accurate study, together with a 
better understanding of lattice effects in the propagators. 
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Figure 3: Cucchieri-Mendes bounds for the lattice data al P — 5.7. 
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